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I.  INTRODUCTION 


The  Air  Force  has  developed  an  Air  Quality  Assessment  Model  (AQAM) 
to  predict  the  dispersion  of  pollutants  emitted  from  typical  Air  Base 
related  sources  (Ref.  1,  2  and  3).  That  model  has  also  been  modified  to 
better  simulate  operations  from  a  Naval  Air  Station  (Ref.  4  and  5).  The 
full  model  consists  of  three  computer  programs ;  Source  Inventory,  Short 
Term  and  Long  Term  Dispersion.  In  the  Short  Term  Model  hourly  averaged 
concentrations  are  predicted  over  a  grid  of  receptors.  Air-base,  air¬ 
craft,  and  off  air-base  (environ)  sources  are  considered  and  are  modeled 
as  point,  area,  or  line  sources  as  appropriate. 

In  the  Short  Term  Model,  the  jet  exhausts  from  aircraft  during  taxi 
and  takeoff  are  treated  as  line  sources  with  no  plume  rise.  Initial 
plume  dimensions  must  be  estimated  as  well  as  the  distance  the  jet 
"penetrates"  the  atmosphere  before  coming  essentially  to  rest  relative  to 
the  ambient  wind.  Stationary  sources  such  as  engine  run  15)  test  stands 
and  test  cells  are  treated  as  points  sources.  Likewise,  the  initial 
plume  dimensions  must  be  estimated  and  a  determination  whether  plume  rise 
is  a  factor  must  be  made.  The  model  then  uses  these  input  data  to  deter¬ 
mine  an  initial  dispersion  coefficient  from  which  a  pseudo  upwind  point 
or  line  source  is  located.  The  sigma's,  or  dispersion  coefficients,  are 
based  on  the  stability  class  of  the  atmosphere  and  the  downwind  distance 
or  travel  time  for  an  elevated,  relatively  low  velocity  source.  There  is 
a  question  then,  as  to  the  applicability  of  these  dispersion  coefficients 
to  turbojet  exhausts  which  are  often  sonic,  have  high  shear  stress/turbu¬ 
lent  mixing  characteristics  and  have  different  turbulence  scales  than 
typical  exhaust  stacks.  Evidence  that  there  is  a  discrepancy  has  been 
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indicated  by  an  Air  Force  study  (Ref.  6)  which  presented  a  pictorial  in¬ 
vestigation  of  the  effects  of  atmospheric  stability  on  a  jet  exhaust. 

It  is  evident  that  qualitative  and  quantitative  data  are  needed  to  char¬ 
acterize  the  dispersion  of  a  turbojet  exhaust  during  ground  operations. 
This  investigation  attempted  to  determine  turbojet  exhaust  dispersion 
as  a  function  of  jet  characteristics  and  wind  direction  under  laboratory 
simulated,  neutrally  stable  atmospheric  conditions. 

The  first  part  of  the  study  involved  the  correct  simulation  of  the 
lower  portion  of  the  atmoshpere  in  an  aeronautical  type  wind  tunnel. 
Published  studies  (Refs.  7,  8  and  9)  have  indicated  that  it  is  not 
sufficient  to  simulate  only  the  mean  velocity  profile  of  the  atmosphere. 
The  turbulent  structure  of  the  atmosphere  must  also  be  correctly  modeled 
simultaneously  with  the  mean  velocity  profile.  Proper  simulation  can  be 
accomplished  following  the  work  of  Sundaram  and  Ludwig  (Ref.  7) ,  tripping 
the  boundary  layer  with  a  suitable  fence  at  the  entrance  of  the  wind 
tunnel  test  section  and  then  letting  the  turbulent  flow  develop  over  a 
length  of  roughness  elements.  Suitable  matching  between  the  boundary 
layer  trip  and  the  roughness  elements  then  results  in  an  area  of  constant 
shear  stress  which  adequately  models  the  atmospheric  flow. 

Once  an  atmoshperic  surface  layer  has  been  suitably  simulated,  a 
turbojet  exhaust  must  then  be  simulated.  If  it  is  assumed  that  heat  and 
mass  have  the  same  turbulent  diffusion  characteristics,  then  the  disper¬ 
sion  rates  of  pollutants  can  be  determined  by  measuring  the  temperature 
distributions  in  the  vertical  and  horizontal  directions. 

In  this  investigation  data  were  obtained  for  varying  angles  of  inci¬ 
dence  to  the  oncoming  wind  for  choked  and  unchoked  nozzle  flow  conditions, 


for  two  jet  temperatures,  and  for  two  conditions  of  nozzle  surface  block¬ 
age.  The  data  were  compared  to  the  dispersion  routines  used  in  the  AQAM 
computer  program. 


ATMOSPHERIC  SIMULATION 


In  order  to  make  meaningful  wind  tunnel  measurements  to  determine 
jet  aircraft  exhaust  dispersion  rates  during  taxi,  idle  and  takeoff 
modes  of  operation  it  was  necessary  to  correctly  model  the  lower  struc¬ 
ture  of  the  atmoshpere.  Sundaram  and  Ludwig  (Ref.  7)  have  published  a 
theoretical  study  of  the  requirements  for  the  modeled  flow  and  an  ex¬ 
perimental  study  of  the  flew  generated  by  wind  tunnel  techniques.  This 
reference  was  used  as  the  principal  guide  in  modeling  the  atmoshperic 
surface  layer  in  a  wind  tunnel  at  the  Naval  Postgraduate  School.  Sup¬ 
plementary  discussions  are  presented  in  Refs.  10-12. 

A.  THEORETICAL  REQUIREMENTS  FOR  ATMOSHPERIC  MODELING 

Reference  13  describes  the  atmoshperic  boundary  layer  as  a  turbulent 
layer  which  is  influenced  by  a  combination  of  actions  due  to  surface 
friction,  thermal  stratification,  and  Coriolis  force.  The  atmospheric 
boundary  layer  is  divided  into  two  distinct  layers.  The  lower  region 
called  the  surface  boundary  layer  extends  up  to  a  nominal  height  of  100 
meters.  The  region  above  the  surface  layer  is  called  the  Ekman  layer 
and  extends  up  to  a  height  of  about  1000  meters.  In  the  surface  layer, 
the  wind  direction  is  generally  considered  to  remain  constant  and  the 
Coriolis  force  is  not  an  important  parameter.  In  an  ideal  environment 
the  near- surf ace  region  is  horizontally  homogeneous  and  free  of  pressure 
gradients.  These  conditions  imply  that  the  shear  stress  can  be  considered 
to  be  independent  of  height  in  the  region  close  to  the  surface.  Surface 
layer  theories  (Ref.  14)  often  employ  the  assumption  that  the  flow  char¬ 
acteristics  close  to  the  ground  can  be  expressed  entirely  in  terms  of 
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conditions  at  the  surface.  Most  atmoshperic  boundary  layers  are  observed 

to  have  no  laminar  sublayer  near  the  surface  and  are  described  to  be  fully 

aerodynamically  rough.  A  turbulent,  neutrally  stable  flow  structure  in 

the  atmospheric  surface  boundary  layer  can  be  described  entirely  by  the 

following  parameters;  the  kinematic  viscosity  v  ,  the  friction  velocity 

uA  ,  and  the  roughness  parameter  z  (Ref.  7).  u.  is  defined  as  being 

o  * 

L 

equal  to  (  t  /p  )  and  z  is  a  length-scale  which  describes  the  irflu- 
o  o 

ence  of  a  rough  surface  on  a  particuJ  ar  flow. 

Using  the  above  parameters,  an  aerodynamically  rough  flow  exists  if 
(Ref.  7) 

ZoU*/v>3*  (1) 

If  zq  is  «  than  the  vertical  dimension  of  the  boundary  layer,  then  the 
mean-velocity  profile  of  the  wind  is  given  in  logarithmic  form  by  the  re¬ 
lationship  (Ref.  7) 

u/u4  -  (1/tc)  In  (z/z  ),  (2) 

o 

where  u  is  the  average  velocity,  z  is  the  vertical  distance  from  the 
surface  and  k  is  Von  Karman's  constant.  This  equation  is  a  form  of  the 
well  known  law  of  the  wall  and  is  indicative  of  the  interaction  between 
the  air  flow  above  a  given  surface  and  that  surface. 

To  accurately  model  the  atmosphere  the  requirements  of  aerodynami¬ 
cally  rough  flow,  horizontal  homogeneity  and  the  absence  of  pressure 
gradients  must  be  adhered  to.  Strict  horizontal  homogeneity,  however,  is 
extremely  difficult  to  obtain  in  the  laboratory.  Partial  simulation, 
therefore,  generally  has  to  be  imposed.  A  complete  discription  of  the 
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effects  of  partial  simulation  is  given  in  reference  7  and  will  not  be 
repeated  here.  The  conclusion  reached  from  reference  7  regarding  partial 
simulation  was  that  a  laboratory  flow  which  was  fully  aerodynamically  rough , 
nearly  horizontally  homogeneous  and  relatively  free  from  pressure 
gradients  was  a  reasonable  model  for  a  neutrally  stable  atmosphere.  The 
generation  of  such  a  flow  is  described  in  the  next  section  of  this  report. 

B.  EXPERIMENTAL  REQUIREMENTS  TO  SIMULATE  ATMOSPHERIC  FLOWS 

Aircraft  jet  exhaust  dispersion  in  the  "near  ground"  environment  was 
the  primary  concern  of  this  investigation;  therefore  only  the  surface  lay¬ 
er  of  the  atmosphere  was  required  to  be  simulated.  As  discussed  previous¬ 
ly,  this  was  possible  since  both  the  atmospheric  flow  and  the  simulated 
flow  are  described  completely  in  terras  of  conditions  at  the  surface  (ie. 
the  parameters  ut  ,  z^  ,  and  v  )  .  The  similarity  requirements  were, 
therefore,  not  a  function  of  the  geostrophic  wind  in  the  real  flow  nor  a 
function  of  the  free  stream  velocity  in  the  simulated  flow.  Horizontal 
homogeneity  cannot  be  strictly  satisfied  since  for  an  experimental  flow 
over  a  flat  plate,  equilibrium  conditions  change  continually  as  a  function 
of  downstream  distance  (Ref.  7)  . 

Since  the  atmosphere  is  assuned  to  exhibit  a  logarithmic  mean  velocity 
profile  and  to  be  horizontally  homogeneous,  this  part  of  the  similarity 
requirement  can  be  verified  by  measuring  the  mean  velocity  profile  at 
several  locations  in  the  wind  tunnel  flow.  In  reference  7  the  turbulence 
intensity  and  the  integral  scale  of  turbulence  were  also  measured  and  com¬ 
pared  to  available  atmospheric  data.  Turbulence  intensity  measurements 
were  also  made  in  this  experiment  and  compared  to  the  results  of  reference 
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7.  Integral  scales  of  turbulence  were  not  measured.  For  virtually 
identical  experimental  set-ups  it  was  assisted  that  if  all  other  experimental 
measurements  compared  favorable  to  those  of  reference  7,  then  the  scales 
of  turbulence  would  also  compare  favorably. 

C.  EXPERIMENTAL  APPARATUS  AND  MEASUREMENTS 

1.  Wind  Tunnel  Modifications 

The  experiment  was  conducted  in  a  low  speed  wind  tunnel  at  the 
Naval  Postgraduate  School  (NPS) ,  Monterey,  CA.  The  wind  tunnel  draws  am¬ 
bient  air  from  the  surrounding  area  into  a  5X5  foot  test  section  approx¬ 
imately  21  feet  long.  In  order  to  be  comparable  with  the  wind  tunnel  used 
in  reference  7,  the  wind  tunnel  at  the  Naval  Postgraduate  School  was  modi¬ 
fied  in  the  following  manner.  A  false  ceiling  was  installed  to  maintain 
an  approximate  zero  pressure  gradient  in  the  axial  direction.  A  two  inch 
high  boundary  layer  trip  constructed  of  wood  was  placed  across  the  entrance 
to  the  test  section.  The  trip  was  followed  by  2.5  X  2.5  feet  square 
sections  of  roughness  elements  which  covered  the  entire  floor  of  the  test 
section.  The  roughness  elements  were  4  sided  regular  pyramids  which  were 
0.75  inches  square  at  the  base  and  0.75  inches  in  height.  They  were  con¬ 
structed  out  of  1.5  inch  thick  fiberboard.  The  tunnel  and  roughness  ele¬ 
ments  are  shown  in  figure  1. 

2.  Pressure  Measurements 

Six,  1/8  inch  diameter  pressure  taps  spaced  4  feet  apart,  were 
mounted  flush  along  the  centerline  of  the  wind  tunnel  wall.  The  differen¬ 
tial  area  required  to  maintain  a  zero  pressure  gradient  in  the  axial  direc¬ 
tion  was  determined  from  an  initial  set  of  pressure  readings.  The  false 
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ceiling  which  implemented  the  required  area  change  was  then  installed  and 


the  pressure  readings  were  again  taken  to  verify  the  required  zero  axial 
pressure  variation. 

3.  Electronic  Equipment 

All  velocity  related  measurements  were  taken  with  a  Thermo  Systems 
Inc.  linearized  hotwire  anemometer  system.  DISA  single  wire  and  cross  win 
probes  and  probe  holders  were  used  with  the  anemometer  system  (fig.  2  and  3) . 
A  DC  digital  voltmeter  and  a  true  RMS  meter  were  used  to  record  the  steady 
and  fluctuating  single  wire  probe  data  respectfully.  A  model  1015C  Thermo 
Systems  Inc.  correlator  was  used  in  conjunction  with  the  cross  wire  outputs 
to  obtain  the  sum  and  difference  of  the  two  signals  (fig.  4) .  The  true 
RMS  meter  was  also  used  to  record  the  cross  wire  outputs  via  the  correlator. 
The  hotwire  probes  were  mounted  under  the  wind  tunnel  such  that  the  probes 
traversed  the  wind  tunnel  in  a  vertical  direction.  A  Spectral  Dynamic's , 
Spectra-Scope  Modal  SD- 330  frequency  counter  was  used  in  conjunction  with 
the  measurement  of  the  centerline  velocity  (fig.  5). 

4.  Wind  Tunnel  Centerline  Velocity  Measurements 

The  maximum  velocity  expected  in  the  wind  tunnel  was  between  40 
and  50  feet  per  second.  Existing  pitot-static  systems  at  NPS  were  not  sen¬ 
sitive  to  this  small  velocity.  Therefore,  the  relationship  between  the 
Strouhal  Number  and  the  Reynolds  Number  was  used  to  determine  the  center- 
line  velocity  (Ref.  15) .  The  Strouhal  number  relates  the  velocity  of  a 
flow  to  the  frequency  at  which  vortices  are  shed  behind  a  circular  cylin¬ 
der  which  is  inmersed  in  that  flow. 


Sn (Strouhal  Number)  =  (freq.  X  diam.  of  cyl.)  /  velocity 


A  0.075  inch  diameter  stainless  steel  wire  was  attached  to  a  collar  which 
was  fitted  over  a  single  hotwire  probe  such  that  the  hotwire  was  in  only 
one  wake  (shed  vortex  sheet)  of  the  cylinder  (fig.  6) .  The  probe,  collar 
and  hotwire  combination  was  placed  perpendicular  to  the  flow  at  the  center- 
line  of  the  tunnel.  With  the  tunnel  on,  the  hotwire  system  was  sensitive  to 
the  velocity  of  the  free  stream  and  the  output  was  observed  and  recorded 
on  the  frequency  counter.  The  observed  frequency  of  shed  vorticies  and 
the  unique  relationship  between  the  Dtrouhal  Number  and  the  Reynolds  Num¬ 
ber  were  used  to  iteratively  determine  the  velocity  of  the  free  stream. 

5.  Velocity  Profile  Measurements 

The  mean  velocity  profile  measurements  were  made  using  the  single 
hotwire  apparatus.  The  hotwire  system  was  calibrated  in  the  linear  mode 
such  that  the  centerline  velocity  was  made  equal  to  one  volt.  The  single 
wire  was  also  calibrated  for  flow  direction  sensitivity.  This  was 
accomplished  by  rotating  the  probe  until  a  maximum  reading  was  observed  on 
the  DC  digital  voltmeter.  It  was  then  assumed  that  the  probe  was  perpen¬ 
dicular  to  the  flow.  All  subsequent  measurements  were  referenced  to  this 
direction.  The  linearized  output  was  averaged  through  a  ten  second  time 
constant  before  being  displayed  on  the  DC  digital  voltmeter.  The  veloci  tv- 
profile  was  obtained  by  recording  the  DC  value  of  the  hotwire  output  as  ii 
traversed  the  wind  tunnel  in  the  vertical  direction.  Vertical  measurements 
were  referenced  to  the  top  of  the  roughness  elements  and  were  measured  to 
within  an  accuracy  of  0.020  inches  with  respect  to  that  reference.  Veloci¬ 
ty  profiles  were  obtained  at  two  axial  locations.  The  locations  were  72 
inches  and  109  inches  downwind  from  the  boundary  layer  trip  at  the  entrance 
to  the  test  section.  The  two  axial  measurements  were  necessary  to  ensure 
the  required  degree  of  uniformity  in  the  jet  dispersion  region. 
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6.  Turbulence/Shear  Stress  Measurements 

Hie  cross  wire  system  described  earlier  was  used  to  measure  the  turbu¬ 
lence  intensities  and  shear  stresses  directly.  The  cross  wire  was  calibra¬ 
ted  such  that  the  linearized  output  for  each  wire  was  of  equal  sensitivity. 
This  was  accomplished  by  first  orienting  one  wire  normal  to  the  flow  and 
then  adjusting  its  sensitivity  to  the  proper  level.  The  probe  was  then 
re-oriented  so  that  the  second  wire  was  normal  to  the  flow  and  its  sensi¬ 
tivity  was  then  matched  to  that  of  the  first  wire.  For  data  acquisition, 
the  probe  was  oriented  axially  such  that  the  longitudinal  turbulence  in¬ 
tensity  was  proportional  to  the  RMS  value  of  the  sum  of  the  cross  wire 
outputs.  Likewise,  for  the  same  orientation,  the  vertical  turbulence  in¬ 
tensity  was  proportional  to  the  RMS  value  of  the  difference  of  the  cross 
wire  outputs.  The  turbulent  shear  stress  is  defined  as  -  pu'w'  .  This 
shear  stress  was  obtained  directly  since  u'w'  is  proportional  to  the 
square  of  the  sum  minus  the  square  of  the  difference  of  the  cross  wire 
outputs.  Like  the  single  wire  data,  the  cross  wire  data  was  recorded  as 
the  probe  was  traversed  in  the  vertical  direction  at  the  same  two  axial 
locations.  For  redundancy,  the  vertical  profile  of  turbulence  intesnity 
was  also  measured  using  the  single  wire  probe. 

D.  RESULTS  AMD  DISCUSSION 

All  the  measurements  obtained  at  the  two  axial  locations,  ie.  72  and 
109  inches  from  the  entrance  of  the  wind  tunnel  test  section,  were  taken 
with  a  wind  tunnel  center  line  velocity  of  approximately  40  feet  per 
second. 
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,  u.  and  z 

1.  _ o 

From  the  law  of  the  wall, 

u/u*  =  5.75  log  z/z  (3) 

o 

Equation  (3}  can  be  rewritten  as 

u/U  =  (uVU  )  5.75  (log  z  -  log  z  ).  (4) 

00  w  oo  Q 

When  u/U^  is  plotted  vs.  log  z, 

slope  =5.75  (u^/u^) .  (5) 

u/U^  vs.  log  z  was  plotted  for  both  axial  positions  and  a  least 
squares  fit  was  obtained  through  the  data  points  (figures  7  and  8) .  The 
slopes  were  measured  and  u4  was  determined  by? 

uA  =  (slope  •  U^/S^S.  (6) 

u^  values  of  2.32  and  2.52  feet  per  second  were  obtained  at  the  72  and 
109  inch  positions  respectively.  Using  the  least  squares  fit  data,  the 
calculated  u*  was  used  to  determine  zq  by; 

log  z  =  log  z  -  u/(5.75  •  u4)  (7) 

o  * 

z  values  of  0.0215  and  0.0326  were  obtained  at  72  and  109  inches 
o 

respectively. 

2.  Vertical  Distribution  of  Shear  Stress 

Since  u*  is  proportional  to  the  shear  stress. 


a  measure  of  the  average  shear  stress  for  each  position  cam  be  obtained 
from  the  mean  velocity  profile  data.  Non-dimensionalizing  by  U^2  ? 

u'w'  /  u«2  *  <u*/Uj2  (9) 


11 


-a 


The  measured  shear  stress,  u'w'  ,  from  the  cross  wire  data,  was  also  non- 

2 

dimenti  Dualized  by  and  plotted  as  a  function  of  the  vertical  dis¬ 

tance  above  the  roughness  elements  for  both  axial  positions.  Figure  9 
presents  the  cross-wire  data.  Also  shown  are  the  values  of  average  stress 
determined  fr  '-.  the  law  of  the  wall. 

Figures  10  and  11  compare  the  results  with  those  presented  in  refer¬ 
ence  7.  It  is  observed  that  the  shear  stress  continued  to  increase  until 
very  near  the  wall  whereas  the  shear  stress  from  reference  7  began  to  de¬ 
cay  further  from  the  wall.  Figure  10  shows  that  the  shear  stress  was 
approximately  constant  for  heights  from  2  to  5  inches  above  the  rough¬ 
ness  elements  at  the  72  inch  position.  At  the  109  inch  position 
(fig.  11) ,  the  shear  stress  decayed  gradually  from  the  near  wall  region. 

In  an  attempt  to  modify  the  near  wall  shear  stress  distribution,  different 
combinations  of  trip  heights  and  wind  speeds  were  tried.  This  proved  to 
be  unsuccessful.  The  mean  velocity  profile  was  observed  to  vary  much  more 
rapidly  in  the  axial  direction  and  the  relatively  constant  shear  stress 
region  obtained  with  the  two  inch  trip  was  destroyed  in  varying  degrees 
depending  on  which  trip  or  what  tunnel  speed  was  used.  There  was  also  a 
somewhat  greater  variation  in  horizontal  homogeneity  in  this  experiment 
relative  to  reference  7.  This  was  reflected  in  the  increased  slope  of 
figure  8  over  that  of  figure  7.  This  was  also  reflected  in  figure  9  by 
noting  that  the  average  shear  stress  differed  by  about  20  percent 
between  the  two  positions.  Moreover,  zq  increased  by  about  50 
percent  over  the  axial  range  compared  to  about  30  percent  for  the  experi¬ 
ment  of  reference  7.  It  was  also  noted  in  reference  7,  however,  that  the 
degree  of  uncertainty  in  the  atmosphere  may  be  much  larger  than  the  30 
percent  figure. 
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3.  Vertical  and  Horizontal  Turbulence  Intensities 

Figures  12  and  13  compare  the  measured  vertical  and  horizontal 
turbulence  intensity  measurements  at  the  two  axial  locations  with  * 
of  reference  7.  In  general,  the  vertical  and  horizontal  data  agre  - 
tremely  well  with  the  data  of  reference  7.  The  horizontal  turbulo  ’  - 
a  little  lower  for  both  positions  in  the  range  of  z/zq  from  1 

40.  This  corresponded  to  a  a  height  of  about  one  inch  above  the-  i . 

elements  and  was  the  area  of  difference  noted  in  figures  10  and  1 1 
Since,  (1)  a  region  of  constant  shear  stress  was  obtained  ir  • 
agreement  with  reference  7,  (2)  the  turbulence  intensities  were  rc-1 

constant  over  the  same  region,  (3)  the  roughness  parameter  z  v - 

o 

parable  to  that  of  reference  7,  and  (4)  equation  (1)  was  satisfied,  •  • 
simulated  atmoshperic  surface  layer  was  considered  to  be  adequate  ■  • 
initial  study  of  jet  exhaust  dispersion. 
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III.  SIMULATED  JET  DISPERSION  IN  A  WIND  TUNNEL 

A.  BACKGROUND 

The  AQAM  model  treats  most  aircraft  emissions  as  finite  line  source 
Aircraft  engine  test  cells  and  run  up  stands  are  treated  as  point  source 
Line  sources  are  modeled  by  a  general  line  source  or  puff  type  model 
(Ref.  1-3) .  Point  sources  are  modeled  by  the  application  of  Gaussi 
plume  theory.  This  theory  represents  the  downwind  concentration  of  ; 
lutants  from  an  elevated  source  as  Gaussian  distributions  in  both  rh-~- 
horizontal  and  vertical  directions.  When  applied  to  aircraft  exhaust..  , 
the  line  source  method  first  assumes  that  the  emissions  from  the  movin'; 
aircraft  come  to  rest  relative  to  the  ambient  wind  after  a  penetration 
length  due  to  the  velocity  effects  of  the  jet.  After  coming  to  rest,  t.-. 
source  of  finite  cross-section  is  generally  segmented.  The  segment?'" 
portions  of  the  source  are  then  assumed  to  behave  as  pseudo-upwind  V.  re¬ 
sources  which  are  dispersed  downwind  in  a  Gaussian  manner.  Point  an j 
line  sources  both  employ  the  same  dispersion  coefficients  and 

which  are  the  standard  deviation  points  of  a  Gaussian  or  normal  di : ' 
tion  curve  and  are  described  by  Turner  in  his  workbook  on  dispersion 
estimations  (Ref.  16).  As  mentioned  earlier,  the  dispersion  coeffic  ■: 
are  based  on  the  atmospheric  stability  class  and  the  time  or  distar.:-.- 
traveled  downwind  from  a  relatively  low  velocity  source.  The  applied  i 
of  these  coefficients  relative  to  a  high  velocity,  horizontally  emitti 
source,  such  as  a  turbojet  engine,  was  investigated  for  various  jet  ;  .. i 
meters. 
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B.  EQUIPMENT 


1.  General  Discussion 

The  investigation  was  conducted  in  the  low  speed  wind  tunnel 
described  in  the  previous  section,  h  bumer/nozzle  assembly  provided  the 
simulated  turbojet  exhaust.  The  temperatures  within  the  exhaust  were 
measured  in  the  vertical  and  horizontal  directions  using  a  thermocouple 
apparatus . 

2.  Bumer/Noyzle 

Oxygen  and  ethylene  were  mixed  and  ignited  at  the  entrance  to  a 
burner  can,  figure  14.  The  hot  gases  were  diluted  with  air  to  provide 
burner  cooling  and  the  desired  flow  rates.  The  burner  exhausted  into  the 
wind  tunnel  through  a  1/8  inch  diameter  stainless  steel  tube.  The  tube 
was  inserted  through  the  tunnel  floor.  It  was  made  with  one  90  degree 
bend  which  aligned  the  exhaust  parallel  to  the  tunnel  floor.  A  nozzle  was 
created  by  reducing  the  exit  diameter  of  the  tubing  to  0.058  inches.  The 
tube  was  positioned  3  inches  above  the  roughness  elements.  This  position 
was  chosen  since  it  was  about  the  center  of  the  area  of  constant  shear 
stress  discussed  in  the  atmospheric  simulation  section  of  this  report. 
Stagnation  pressure  and  temperature  probes  were  used  to  determine  the 
nozzle  exit  conditions.  The  nozzle  entered  the  wind  tunnel  at  a  position 
72  inches  downstream  from  the  boundary  layer  trip  described  in  the  previ¬ 
ous  section. 

3.  Temperature  Measuring 

Horizontal  and  vertical  tenqperature  profiles  were  obtained  using  a 
chromel-alumel  thermocouple  probe.  The  thermocouple  was  referenced  to  an 
Omega  electronic  Miniature  Cold  Junction  (MCJ) .  The  nozzle  and 
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thermocouple  system  is  shown  in  figure  15.  A  probe  holder  and  a  traverse 


system  permitted  positioning  of  the  thermocouple  probe  in  varying  axial, 
transverse,  and  vertical  locations.  Hie  traverse  system,  an  electrical- 
mechanical  device  mounted  under  the  v/ind  tunnel  (fig.  16)  ,  facilitated 
positioning  in  the  vertical  and  transverse  direction.  The  probe  holder, 
mounted  on  the  traverse  system,  permitted  the  thermocouple  probe  to  be 
positioned  in  the  axial  direction.  Two  X-Y  plotters  (fig.  17)  were 
used  to  record  the  temperature  profiles  as  a  function  of  downwind  distance. 
One  plotter  was  used  for  the  vertical  temperature  profile  and  the  other 
for  the  horizontal  temperature  profile. 

C.  DISCUSSION  OF  RESULTS 

1.  Introduction 

The  dispersion  data  from  the  experimental  jet  exhaust  were  used  to 
determine  1)  if  the  plume  dispersion  was  Gaussian,  2)  the  sensitivity 
of  the  jet  dispersion  to  the  wind  direction,  3)  the  sensitivity  of  the 
jet  exhaust  dispersion  to  nozzle  exit  Mach  number  (choked  and  unchoked) 
under  no-wind  and  wind  conditions,  4)  the  sensitivity  of  the  jet  disper¬ 
sion  to  upwind  nozzle  surface  blockage  and  5)  the  effect  of  increased 
jet  exhaust  exit  temperature  on  plume  rise.  The  first  four  items  were  in¬ 
vestigated  with  an  exhaust  stagnation  temperature  of  550  degrees  F  while 
the  stagnation  tenperature  for  the  plume  rise  investigation  was  increased 
to  700  degrees  F.  For  choked  flow  the  exhaust  stagnation  pressure  was  32 
psia.  For  unchoked  flow  it  was  22  psia. 

2.  Attainment  of  Gaussian  Dispersion 

A  Gaussian  or  normal  distribution  with  a  mean  value  of  zero  is 
defined  by  the  equation 
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y  =  { 1/ ( )  >  exp[-.5(x/a) 2] 


(10) 


A  standard  normal  distribution  (fig.  18)  is  a  normal  distribution  in  which 
o  is  set  equal  to  one.  In  this  case,  the  maximum  value  of  the  ordinate 
of  equation  10  is  l/t/JiT  or  0.3989.  The  value  of  the  ordinate  at  one 
standard  deviation,  where  x  =  a  ,  of  equation  10  is  0.242.  To  determine 
whether  the  measured  temperature  profiles  were  Gaussian  the  profiles  were 
standardized  in  the  following  manner.  For  any  measured  temperature  profile 
the  peak  temperature  (Tmax)  was  assumed  to  be  at  the  mean  of  the  distri¬ 
bution.  That  temperature  was  scaled  such  that  Tmax/C  was  equal  to  0.3989. 
The  ordinate  values  of  the  temperature  profile  were  then  scaled  by  the 
factor  C.  The  ordinate  value  of  the  one  a  point,  T^  ,  of  the  profile 
was  determined  by;  Ta  =  0*  0.242.  At  this  temperature  the  width  of  the 
temperature  plume  is  equal  to  2o  .  The  standard  deviation  which  was 
determined  in  this  manner  was  used  to  non-dimensionalize  or  standardize 
the  abscissa  of  the  temperature  profile.  The  standard  deviations  are  ta¬ 
bulated  in  Tables  I,  II  and  III  for  each  of  the  conditions  of  the  experi¬ 
ment.  The  standardized  data  points  were  then  plotted  against  a  standard 
normal  distribution  curve  and  compared  for  closeness  of  fit.  Figures  19 
through  24  compare  the  vertical  and  horizontal  temperature  profiles  for  a 
choked  nozzle  flow  to  the  standard  normal  distribution  curve.  The  jet 
exhaust  was  oriented  at  zero  degrees  incidence  to  the  free  stream  wind 
direction.  The  wind  tunnel  centerline  velocity  was  40  feet  per  second. 

The  different  sets  of  data  represent  measurements  taken  at  different  axial 
positions  downwind  from  the  jet  exhaust.  The  downwind  distance  was  ex¬ 
pressed  in  terms  of  jet  diameters.  The  plume  was  near  Gaussian  after 
approximately  eight  jet  diameters  downstream  from  the  exit  plane  of  the 
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jet  for  both  the  horizontal  and  vertical  profiles.  The  distanoe  required 
for  the  profile  to  become  Gaussian  appeared  to  be  the  result  of  the  expan¬ 
sion  process  of  the  fluid  as  it  exited  the  nozzle.  Similar  plots  were  ob¬ 
tained  for  the  other  zero  incidence  conditions  shown  in  Tables  I  and  ll. 

For  each  case  the  plume  exhibited  near  Gaussian  properties  at  about  eight 
jet  diameters  downstream.  The  dashed  line  on  Tables  I  and  II  indicate  the 
position  at  which  the  plume  exhibited  Gaussian  characteristics. 

As  the  nozzle  was  rotated  to  the  wind  direction,  the  jet-wind  inter¬ 
actions  caused  the  horizontal  profile  to  be  negatively  skewed.  The  degree 
of  skewness  diminished  as  a  function  of  downwind  distance  as  the  effect  of 
the  jet  exit  velocity  was  overcome  by  the  force  of  the  wind.  At  the  point 
where  the  wind  was  the  dominant  dispersion  factor  (ie.  past  the  penetration 
length)  the  plume  again  exhibited  Gaussian  characteristics.  Figure  25  com¬ 
pares  the  data  obtained  at  an  axial  distance  of  1.72  jet  diameters  for 
different  angles  of  incidence  to  the  wind.  The  angles  investigated  were 
29,  45,  67  and  90  degrees  and  the  data  were  non-dimensionalized  in  the 
same  manner  as  discussed  previously.  The  degree  of  skewness  in  figure  25 
is  evident  for  each  angle  of  incidence  when  data  are  compared  to  the  stan¬ 
dard  normal  distribution  curve.  Similar  data  for  different  downwind  sta¬ 
tions  are  plotted  on  figures  26  through  32.  The  vertical  temperature 
profiles  exhibted  no  skewness  and  were  very  nearly  Gaussian  at  each  data 
position.  The  a's  which  were  obtained  for  the  vertical  and  horizontal 
dispersions  are  tabulated  in  Table  III.  The  downwind  distance  at  which 
the  plume  became  Gaussian  is  indicated  for  each  angle  by  a  dotted  line. 
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3.  Sensitivity  of  Exhaust  Dispersion  to  Angle 
of  Incidence  to  the  Wind 

A  spatial  representation  of  the  jet  exhaust  as  a  function  of  dcwn- 
wind  and  horizontal  distances  from  the  centerline  of  the  nozzle  is  plot¬ 
ted  on  figure  33  for  jet  exhaust  angles  of  incidence  to  the  wind  of  29 , 

45,  67  and  90  degrees.  As  expected,  as  the  angle  of  incidence  was  in¬ 
creased,  the  length  normal  to  the  wind  direction  where  the  dispersion 
tribution  was  essentially  uniform  (ie.  sigma  approaches  infinity)  incrc 
(in  a  nonlinear  relationship,  fig.  34).  The  downwind  distance  at  whi 
the  plume  became  approximately  uniformly  distributed  was  extremely  sensi¬ 
tive  to  the  angle  of  incidence.  The  downwind  distance  decreased  quite 
rapidly  as  the  angle  of  incidence  was  increased  in  an  almost  linear  re¬ 
lationship  (fig.  34).  From  Table  III  it  can  also  be  seen  that  the 
corresponding  sigmas  increased  considerably.  Quantitative  data  was 
attempted  for  angles  of  incidence  greater  than  90  degrees.  The  prob.. 
could  not  be  translated  further  upstream  than  the  90  degree  jet  orienta¬ 
tion.  At  this  position,  the  exhaust  of  the  jet  was  completely  disperse. 
Qualitatively  then,  as  the  angles  of  incidence  were  increased,  the  exit 
jet  momentum  and  wind  action  acted  together  to  rapidly  disperse  the  ex¬ 
haust  in  the  downwind  direction.  The  dispersion  coefficients  were  in¬ 
creased  which  resulted  in  a  decreased  downwind  distance  where  the  plume 
became  uniformly  distributed. 

4.  Sensitivity  of  Dispersion  Rate  to  Nozzle  Exit  Mach  Number 
(choked  and  unchoked) ,  for  Wind  and  No-wind  Conditions 
Table  I  presents  the  results  in  terms  of  the  size  of  the  one 

standard  deviation  point  of  the  plume  in  jet  diameters  as  a  function  of 
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downwind  distance  from  the  nozzle  exit.  It  can  be  seen  that  the  plume  di¬ 
mensions  were  relatively  insensitive  to  varying  exit  conditions  (choked, 
unchoked  and  high  temperature) .  The  data  for  the  unchoked  nozzle-with-wind 
dispersed  about  the  same  as  did  the  choked  flow.  For  choked  flow  the 
Mach  mxnber  at  the  exit  was  1.0.  For  the  unchoked  flow  the  Mach  number 
was  about  0.8.  Therefore,  both  had  high  velocities  relative  to  a  typical 
stack  exhaust. 

The  dispersions  from  a  choked  nozzle  with  and  without  wind  were  also 
measured.  The  two  sets  of  data  for  both  the  horizontal  and  vertical  di¬ 
rections  are  plotted  in  figures  35  and  36.  It  is  seen  that  the  disper¬ 
sion  of  the  plume  became  affected  by  the  wind  structure  at  about  30  jet 
diameters  downwind  from  the  nozzle  exit.  Between  the  nozzle  exit  and 
about  30  jet  diameters  downstream,  the  no-wind  and  wind  data  compare 
fairly  well  with  each  other.  Further  downwind  from  that  position,  however, 
the  two  curves  diverge.  It  was  apparent  then  that  the  exit  conditions  of 
the  nozzle  (ie.  the  penetration  length)  governed  the  dispersion  of  the 
plume  for  about  the  first  30  jet  diameters  downwind,  and  thereafter  the 
dispersion  was  primarily  a  function  of  the  environmental  conditions. 

5.  Effect  of  Nozzle  Surface  Blockage 

To  investigate  the  effect  of  flow  disturbances  (caused  by  the 
tubing  upwind  of  the  nozzle)  on  the  dispersion  of  the  jet  exhaust,  two 
different  non-heat  conducting  disks  were  separately  fitted  on  the  nozzle 
tubing.  Hiese  blockage  devices  increased  the  upstream  area  relative  to 
the  1/8  inch  tubing  by  a  factor  of  15  and  75  respectfully.  Table  II 
contains  the  results  and  it  was  observed  that  upstream  blockage  did  not 
greatly  affect  the  dispersion  rates.  These  upstream  blockage  effects 
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were  apparently  dissipated  by  the  time  that  the  jet  had  penetrated  approx- 

i 

imately  30  jet  diameters  aft  of  the  nozzle.  These  data  indicate  that  air¬ 
craft  configuration  may  not  significantly  affect  the  jet  dispersion  for 
the  slowly  moving  taxi  operation. 

6.  Sensitivity  to  Increased  Exit  Temperature 

As  stated  previously,  the  effect  of  exit  temperature  had  little 
effect  on  the  plume  dimensions  as  a  function  of  downwind  distance.  An 
investigation  of  the  plume  rise  was  conducted  by  measuring  the  vertical 
distance  between  the  peak  point  of  the  temperature  profile  and  the  center- 
line  of  the  nozzle  for  two  exhaust  temperatures.  The  two  sets  of  data  are 
compared  as  a  function  of  downwind  distance  in  figure  37.  Because  of  the 
small  distances  it  was  difficult  to  make  accurate  quantitative  conclusions 
with  respect  to  the  plume  rise.  From  the  data,  however,  it  appeared  that 
the  plume  rise  was  about  one  jet  diameter  per  100  degrees  F  change  in  exit 
temperature.  Equally  as  interesting  as  the  plume  rise  was  what  appeared 
to  be  a  small  degree  of  looping  caused  by  the  turbulent  structure  of  the 
atmosphere.  This  phenomenon  is  shown  in  figure  37  as  a  dotted  line. 

D.  COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  AQAM 

The  line  source  dispersion  enployed  for  aircraft  exhaust  jets  in  AQAM 
first  assumes  that  the  jet  exhaust  has  a  penetration  length  equal  to  140 
meters  (ie.  approximately  140  jet  diameters).  The  penetration  length  is 
that  length  where  the  velocity  effects  of  the  exhaust  cease  to  be  a  fac¬ 
tor  in  dispersing  the  exhaust.  Thus,  the  exhaust  essentially  comes  to  rest 
with  respect  to  the  ambient  air  mass.  At  this  distance  then,  the  exhaust 
is  being  dispersed  entirely  by  the  action  of  the  wind. 

In  AQAM,  estimates  of  the  lateral  and  vertical  dimensions  of  the 
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aircraft  plane  are  made  at  the  penetration  length.  The  lateral  width  w  , 

and  vertical  height  Az  then  define  the  dimensions  of  a  line  source  of 

length  L  .  The  lateral  width  is  assumed  to  remain  constant  over  the 

length  of  the  line.  Once  the  line  source  has  been  defined,  the  line  is 

generally  segmented  and  dispersed  downwind  using  the  Gaussian  theory  of 

dispersion.  That  is;  from  the  dimensions  of  the  line  source  an  initial 

and  Ozo  are  determined  from  empirical  relationships  =  w/2.4  and 

o  =  A*/2.4  respectively.  With  these  values  of  a  and  a  an  up- 
zo  yo  zo 

wind  distance  for  a  pseudo  upwind  line  source  is  determined  using  the 

a  and  o  vs.  distance  or  time  charts.  The  pollution  concentrations 
yo  zo 

are  then  dispersed  from  the  determined  line  source  in  a  Gaussian  manner. 

If  w  and  Az  are  not  specified  inputs,  the  AQAM  model  uses  default 
values  of  20  meters  and  8  meters  respectfully. 

The  penetration  length  determined  in  this  investigation  for  zero 
degrees  orientation  to  the  wind  was  approximately  30  jet  diameters,  con¬ 
siderably  shorter  than  the  140  jet  diameters  used  in  AQAM.  The  30  jet 
diameter  figure  was  further  investigated  using  the  dispersion  data  from 
tests  with  the  nozzle  oriented  at  other  than  zero  degrees  to  the  wind 
(fig.  38).  The  30  jet  diameter  penetration  distance  determined  from  the 
axial  conditions,  indicated  by  an  arc  in  figure  38,  was  compared  to  the 
point  where  the  dispersion  curves  exhibited  no  skewness  (ie.  where  the 
distributions  were  Gaussian  normal  to  the  wind  direction) .  It  was 
assumed  that  beyond  this  position  the  jet  effects  could  be  neglected  and 
that  the  jet  would  be  disperesed  by  the  ambient  wind  only.  This 
assumption  is  not  entirely  correct  since  the  plume  continued  to  exhibit 
horizontal  displacements  normal  to  the  wind.  Along  the  jet  axis  the 
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penetration  length  increased  as  a  function  of  incidence  to  the  wind  from 
about  28  to  58  jet  diameters,  still  considerably  less  than  the  140  jet  di¬ 
ameters  assumed  by  AQAM.  This  relationship  is  plotted  in  figure  39  and 
appears  to  be  approximately  linear  between  29  and  90  degrees  relative  to 
the  axial  orientation.  The  distance  downwind  from  the  nozzle  where  the 
plume  became  Gaussian,  however,  remained  approximately  constant  at  25-30 
jet  diameters.  It  was  also  noted  that  the  corresponding  o's  at  the 
Gaussian  point  for  each  angle  became  wider  with  increased  angles  of  inci¬ 
dence  due  to  more  turbulent  mixing.  This  relationship  is  plotted  in 
figure  40.  With  the  limited  number  of  data  points  available  it  appears 

that  the  relationship  follows  a  smooth  curve,  a  and  o  and  the  jet 

y  z 

penetration  length  seem  to  be  predictable  between  zero  and  90  degrees  of 
incidence  to  the  wind  for  neutral  stability  and  these  nozzle  exit  con¬ 
ditions.  These  values  could  be  used  in  AQAM  to  determine  w  and  Az 
instead  of  estimating  them  or  using  the  default  values.  The  relationship 
shown  in  figure  39  could  also  be  used  as  input  for  jet  penetration  length. 
However,  since  the  penetration  lengths  are  small,  it  may  be  better  to 
neglect  this  effect  entirely. 

A  comparison  of  the  rates  of  dispersion  for  zero  degrees  incidence 

to  the  wind  (i.e.  the  slopes  of  the  o  and  c  vs.  distance  curves)  was 

y  z 

made  between  those  for  a  neutrally  stable  atmoshpere  (class  D)  in 
reference  1  and  the  experimental  data.  As  can  be  seen  from  figures  35 
and  36  the  ejq>erimental  dispersion  rates  were  greater  than  the  class  D 
rate  used  by  reference  1  for  a  neutrally  stable  atmosphere.  Reference  1 
indicated  that  this  may  be  expected.  The  dispersion  rates  were  obtained 
from  data  compiled  from  typical  low  velocity  stack  emissions  and  not  from 
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high  velocity  sources.  In  general  terms,  it  appears  that  for  the  condi¬ 
tions  of  this  experiment,  a  more  realistic  representation  of  the  actual 
dispersion  rate  could  be  obtained  if  the  atmospheric  stability  class  was 
decreased  by  one. 

For  a  vertical  stack,  AQAM  also  estimates  the  plume  rise  due  to 
thermal  effects  by  one  of  two  relationships.  These  relationships  were 
investigated  to  see  if  they  could  be  used  to  predict  the  plume  rise  ob¬ 
served  in  the  experiment.  Both  Holland  and  Carson-Moses  formulas 
were  found  to  predict  insignificant  plume  rise  for  the  conditions  of  this 
experiment. 
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IV.  SUMMARY  AND  RECOMMENDATIONS 

in  summary,  a  neutrally  stable  atmospheric  surface  layer  was  simu¬ 
lated  in  a  low  speed  wind  tunnel.  Horizontal  and  vertical  temperature 
profiles  aft  of  a  converging  nozzle  were  measured,  lhe  dispersion  rates 
from  the  simulated  jet  exhaust  were  measured  for  different  nozzle 
conditions  and  compared  to  the  dispersion  rates  used  in  the  AQAM  Model. 

Shorter  jet  penetration  lengths  and  more  rapid  jet  dispersions  were 
observed.  Initial  plume  dimensions  were  found  to  vary  significantly  with 
jet  orientation  to  the  ambient  wind  direction  and  some  plume  rise  was 
observed. 

It  is  recommended  that  additional  studies  be  made  for  conditions  of 
other  than  neutral  stability.  Measurements  made  in  less  stable  surface 
layers  would  provide  a  better  understanding  of  the  plume  rise  and  looping 
phenomena. 
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3.44 
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4.59 
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4.74 

4.90 

5.17 

86.2 

6.61 
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7.32 

7.50 

7.75 
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9.19 

9.19 

9.19 

11.20 

11.30 
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Table  II 


1-Sigoa  Plume  Width  of  the  Vertical  and  Horizontal 
Temperature  Profiles  at  Zero  Degrees 
Incidence  to  the  Wind  for  Different 
Nozzle  Blockages 
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Figure  1 

Low  Speed  Wind  Tunnel  and  Roughness  Elements 
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Figure  8 

Velocity  Profile  at  109  inches  Downwind 
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Figure  11 

Vertical  Distribution  of  Shear  Stress  at 
109  Inches  Downwind  Compared  to  the  Data  of  Ref.  7 


Figure  12 

Comparison  of  Turbulence  Intensities  at 
2  inches  Downwind  to  the  Data  of  Ref.  7 


Comparison  of  Turbulence  Intensities  at 
109  inches  to  the  Data  of  Ref.  7 


Figure  14 
Burner  Can 
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Figure  15 

Exhaust  Nozzle  and  Thermocouple  Probe 


Figure  lb 

Thermocouple  Traverse  System 
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Comparison  of  Vertical  Distribution  of 
Temperature  to  a  Standard  Normal  Distribution 
for  Zero  Degrees  Angle  of  Incidence 
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Figure  22 

Comparison  of  Horizontal  Distribution  of 
imperature  to  a  Standard  Normal  Distribution 
lor  Zero  Angle  of  Incidence 
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Comparison  of  Horizontal  Distribution  of 
Temperature  to  a  Standard  Normal  Distribution 
for  Zero  Angle  of  Incidence 
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Figure  25 

Comparison  of  Horizontal  Distribution  of  Temperature  to  i 
Standard  Normal  Distribution  for  Angles  of  Incidence  at 
1.72  Jet  Diameters  Downwind  from  the  Nozzle 


of  Incidence 


Figure  27 

Comparison  of  Horizontal  Distribution  of  Temperature  to  ■ 
Standard  Normal  Distribution  for  Angles  of  Incidence  at 
17.2  Jet  Diameters  Downwind  from  the  Nozzle 


Comparison  of  Horizontal  Distribution  of  Temperature  to 
Standard  Normal  Distribution  for  Angles  of  Incidence  a 
86.2  Jet  Diameters  Downwind  from  the  Nozzle 
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Figure  34 

Distance  Downwind  and  Normal  to  Wind  from  the  Nozzle 
Where  the  Plume  is  Uniformly  Distributed 
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Vertical  Plume  Width  vs.  Downwind  Distance 
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Figure  37 

e  Rise  vs.  Downwind  Distance 


Jet  Penetration  vs.  Downwind  and  Horizontal  Distance  from 
the  Nozzle  for  Different  Angles  of  Incidence  to  the  Wind 


Jet  Axis  Penetration  Distance  (JD) 
Figure  39 


Jet  Penetration  Distance  along  the  Jet 
Axis  vs.  Angle  of  Incidence 
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Figure  40 

1-Sigma  Width  of  Plume  At  Gaussian  Position 
vs.  Angle  of  Incidence  to  the  Wind 
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